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Abstract 

We derive a multiconfigurational time-dependent Hartree theory for systems with particle con- 
version. In such systems particles of one kind can convert to another kind and the total number 
of particles varies in time. The theory thus extends the scope of the available and successful mul- 
ticonfigurational time-dependent Hartree methods - which were solely formulated for and applied 
to systems with a fixed number of particles - to new physical systems and problems. As a guiding 
example we treat explicitly a system where bosonic atoms can combine to form bosonic molecules 
and vise versa. In the theory for particle conversion, the time-dependent many-particle wavefunc- 
tion is written as a sum of configurations made of a different number of particles, and assembled 
from sets of atomic and molecular orbitals. Both the expansion coefficients and the orbitals form- 
ing the configurations are time-dependent quantities that are fully determined according to the 
Dirac-Prenkel time-dependent variational principle. By employing the Lagrangian formulation of 
the Dirac-Prenkel variational principle we arrive at two sets of coupled equations of motion, one for 
the atomic and molecular orbitals and one for the expansion coefficients. The first set is comprised 
of first-order differential equations in time and nonlinear in-general integrodifferential equations in 
position space, whereas the second set consists of first-order differential equations with coefficients 
forming a time-dependent Hermitian matrix. Particular attention is paid to the reduced density 
matrices of the many-particle wavcfunction that appear in the theory and enter the equations of 
motion. There are two kinds of reduced density matrices: particle-conserving reduced density 
matrices which directly only couple configurations with the same number of atoms and molecules, 
and particle non-conserving reduced density matrices which couple configurations with a different 
number of atoms and molecules. Closed-form and compact equations of motion are derived for 
contact as well as general two-body interactions, and their properties are analyzed and discussed. 

PACS numbers: 31.15.xv, 05.30.Jp, 05.30.Fk, 03.65.-w 
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I. INTRODUCTION 



The exploration of quantum dynamics of many-particle systems is a fundamental and 
on-going challenge of many branches in physics [1-6]. The equation of motion governing the 
evolution of quantum particles is, in many cases, the well-known time-dependent Schrodinger 
equation. Solving the time-dependent Schrodinger equation for many-particle systems can 
rarely be made analytically or exactly, which renders efficient approximations a must. 

The multiconfigurational time-dependent Hartree method (MCTDH) [7, 8], which has 
been developed in the past two decades, is considered at present the most efficient wave- 
packet propagation approach [9] and has successfully and routinely been used for multi- 
dimensional dynamical systems consisting of distinguishable degrees-of-freedom, such as 
molecular vibrations, see Ref. [10-19]. The main idea behind the MCTDH method is to 
expand the time- dependent many-body wavefunction of distinguishable particles by time- 
dependent configurations that are assembled from time- dependent orbitals (one-body func- 
tions) and optimized according to the Dirac-Prenkel variational principle [20, 21]. In this 
way, a much larger effective subspace of the many-particle Hilbert space can be spanned 
in practice in comparison to multiconfigurational expansions with stationary configurations. 
By grouping several "elementary" degrees-of-freedom together and treating them as "gen- 
eralized" particles, the efficiency of the MCTDH algorithm increases [10, 11]. Choosing to 
use MCTDH itself to propagate multi-dimensional "generalized" particles has led to the 
idea of cascading [13]. Finally, expanding the time-dependent orbitals themselves by other 
time-dependent orbitals, and so on, putting the resulting time-dependent expansion under 
the Dirac-Prenkel variational principle, leads to the multi-layer formulation of the MCTDH 
theory [22] which further increases the efficiency of the MCTDH method for larger, complex 
systems. The MCTDH can be applied to systems of identical particles. In this direction, 
we would like to mention that the MCTDH approach has very successfully been employed 
to unveiling fundamental physics of few-boson systems [23-28] on the numerically-exact 
many-body level. 

A new branch of MCTDH-based methods has emerged after it had been reahzed that, 
to effectively treat the dynamics of more than a handful identical particles, it is essential to 
use their quantum statistics, Fermi-Dirac or Bose-Einstcin, to eliminate the large amount of 
redundancies of coefficients in the distinguishable-particle multiconfigurational expansion of 
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the MCTDH wavefunction. First, taking explicitly the antisymmetry of the many-fermion 
wavefunction to permutations of any two particles into account, the fermionic version of 
MCTDH - MCTDHF - was independently developed by several groups [29-31]. Shortly 
after, the bosonic version of MCTDH - MCTDHB - was developed in [32, 33]. This ad- 
vancement is in particular valuable since very-many bosons can reside in only a small number 
of orbitals owing to Bose- Einstein statistics, thereby allowing the successful and quantitative 
attack of the dynamics of a much larger number of bosons with the MCTDHB theory. For 
applications of MCTDHF to the many-body dynamics of at-present few-fermion systems 
with or without external laser field see Refs. [34-40] , and for applications of MCTDHB for 
the many-body dynamics of repulsive and attractive bosonic systems Rcfs. [32, 41-43]. We 
mention that Ref. [43] has combined optimal control theory with MCTDHB. 

Five decades ago, in his seminal paper, Lowdin defined the reduced density matrices of 
many-fermion wavef unctions [44]. Since then, reduced density matrices and, in particular, 
reduced two-body density matrices is a fruitful and vivid research area including theory and 
applications in electronic structure of molecules, quantum phase-transitions, and ground- 
state nuclear motion [45-51]. In the present context, reduced one- and two-body density 
matrices were first used to derive the stationary many-body states within the general vari- 
ational theory with complete self-consistency for trapped bosonic systems - the multicon- 
figurational Hartree for bosons (MCHB) [52]. Later on, the MCTDHF and MCTDHB were 
formulated in a unified manner, making use of the reduced one- and two-body density ma- 
trices of the time-dependent many-body wavefunction [53]. Finally, treating mixtures of 
two kinds of identical particles in a unified manner, and utilizing the reduced one- and two- 
body density matrices of the mixture's wavefunction, a multiconfigurational time-dependent 
Hartree method for Fermi- Fermi (MCTDH-FF), Bose-Bose (MCTDH-BB) and Bose- Fermi 
(MCTDH-BF) mixtures has been derived [54]. 

The multiconfigurational time-dependent Hartree method and its versions specified for 
identical particles and mixtures are particle- conserving many-body propagation theories. 
Namely, they were solely formulated for and applied to systems with a fixed number of 
particles. Conceptually, they aim at describing systems of coupled degrees-of-freedom or 
interacting particles which have first-quantization Hamiltonian. This brings us to the theme 
of the present work, which is to derive a multiconfigurational time-dependent Hartree the- 
ory for systems with particle conversion. In such systems particles of one kind can convert 
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to another kind and the total number of particles varies in time. Hence, they are gen- 
erally represented by a phenomenological second-quantized Hamiltonian which includes a 
conversion term. Doing so, we extend the scope of the available and successful multicon- 
figurational time- dependent Hartree method and its versions specified for identical particles 
and mixtures to new physical systems and problems. We abbreviate the multiconfigurational 
time-dependent Hartree theory for systems with particle conversion by MCTDH- conversion. 

As a concrete and guiding example for a many-body system with particle conversion 
and without loss of generality, we consider explicitly the conversion of bosonic atoms (a) to 
bosonic molecules (m) via the 'reaction' 2a ^ m, which has been a system of tremendous 
theoretical and experimental interest in quantum-gas physics [55-85]. 

An effective quantum-ficld-theory-based Hamiltonian for atomic and molecular Bose- 
Einstein condensates (BECs) coupled by conversion was first put forward by Drummond et 
al. in [55]. In [56], a proposition that a molecular BEC could be produced by coherent pho- 
toassociation was made and a phenomenological two-mode Hamiltonian to describe this pro- 
cess was suggested. A microscopic theory to derive the many-body Hamiltonian of bosonic 
atoms and molecules with conversion and the respective Gross- Pit aevskii theory with con- 
version were put forward in [57], also sec [59]. In [58], a coupled system of Gross-Pitacvskii 
equations with conversion and deactivation-rate (dissipation) terms has been derived. The 
validity of the two-mode approach for conversion, at least in the homogeneous system, was 
questioned in [62], where dissociation of molecules to other than the ground atomic mode 
signifies that one needs to go beyond the two-mode approximation. The importance of pair 
correlations in the dynamics of resonantly-coupled atomic and molecular BECs, leading to 
significant deviations from the respective Gross-Pitaevskii theory, was put forward in [63]. 
That even in the perfect two-mode limit the mean-field theory with conversion can fail, 
because of strong particle-particle entanglement near the dynamically unstable molecular 
mode, was reported in [65] . A proposition to create a molecule BEC from an atomic Mott- 
insulator phase with exactly two bosons per lattice site was made in [68] . A full microscopic 
theory to derive the Hamiltonian of atoms and molecules with the conversion term from 
the microscopic particle-conserving Hamiltonian of a homogeneous gas of identical bosonic 
atoms with two internal states was given in [69], also see [74]. Quantum phase transitions 
and effects of rotations in homogeneous systems of atomic and molecular BECs with con- 
version have been discovered in [75, 77] and [76], respectively. Finally, in a (harmonic) trap. 
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confinement effects on the stimulated dissociation (effective conversion rate) of a molecular 
to an atomic BEC were recently found in [83] , and unique phases (vortex configurations) of 
rotating interacting atomic- molecular BECs in [85]. 

Molecules were first produced from and identified in a ^^Rb BEC by Wynar et al. [60]. 
Soon after, photoassociation of ultracold sodium molecules in an atomic BEC was made 
[66] . Atomic-molecular coherence in a BEC [made of ®^Rb atoms] was first achieved in [67] . 
A pure molecular quantum gas produced from an atomic cesium BEC was reported in [71] , 
and a quantum-degenerate gas of sodium molecules in [72] . More recently, with ^'''Rb atoms 
in the Mott- phase of optical lattices, state-selective conversion of atoms to molecules [78], 
following the theoretical proposition in [68], and atom- molecule Rabi oscillations [84] have 
been observed. 

Finally, systems with particle conversion can involve of course fermions, in the cold- 
atom world - see the reviews [82, 86] and references therein - and beyond it. In the latter 
context, it is gratifying to mention the Priedberg-Lee model of superconductivity, describing 
the conversion of two electrons to a single Cooper-pair and vice versa by a boson-fermion 
Hamiltonian with a phenomenological conversion term [87, 88]. 

Let us return now to MCTDH- conwerszon, and put it in the particular context of in- 
teracting atomic and molecular BECs with conversion. For the explicit scenario of the 
conversion 'reaction' 2a^m dealt with throughout this work, the theory shall be referred to 
as MCTDH-[2a^m]. MCTDH-[2a^m], as its particle-conserving predecessors [7, 8, 12, 29- 
33, 53, 54], is intended for systems with a finite number of interacting particles, typically 
trapped in an external potential. As a first step, we extend or "merge" two theoretical 
approaches much in use in the literature: the Gross-Pitaevskii theory with conversion and 
the two-mode approximation, see, e.g., Refs. [57, 59, 61, 64, 81, 85] and Refs. [56, 65, 80], 
respectively. This results in a fully variational theory where the two modes - the atomic and 
molecular orbitals - and each and every expansion coefficient in the two-mode many-body 
wavefunction are fully optimized - the orbitals in time and space and the expansion coeffi- 
cients in time - according to the Dirac-Frenkel variational principle [20, 21]. Our main aim is 
to go beyond any two-mode description of the atomic-molecular coupled system and present 
a fully- variational multiconfigurational time-dependent many-body theory for bosonic atoms 
and molecules coupled by conversion - the MCTDH-[2a^m] theory. 

The structure of the paper is as follows. We open in section II with the many-body 
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Hamiltonian of the system of atoms and molecules with conversion. In section III we consider 
as mentioned above a specific case of interest, the fully- variational theory where there are 
one atomic and one molecular orbitals. This specific theory will be referred to as conversion 
mean field. Next, section IV is devoted to the general theory. Both time-dependent as well as 
time-independent theories are presented. Finally, in section V wc put forward a summary and 
concluding remarks. Complementary derivations and relevant matrix elements are deferred 
to and collected in appendices A and B, respectively. 



II. THE MANY-BODY HAMILTONIAN OF INTERACTING ATOMS AND 
MOLECULES WITH CONVERSION 



As a concrete example for a many-body system with particle conversion and without loss 
of generality, we consider a system of bosons which will be referred to as atoms (a) and 
their conversion to another type of bosons which will be referred to as molecules (m) via 
the 'reaction' 2a ^ m. The many-body Hamiltonian of the coupled atom-molecule system 
is taken from the literature of cold-atom physics [55, 57, 69, 76] and is written for our needs 
as a sum of four terms: 



(1) 



The first three terms are particle-conserving terms and together describe a mixture of two 
kinds of interacting bosonic particles; a (atoms) and m (molecules): 



dr 



*l(r)/.('^)(r)*„(r) + i f rfr'*t (r)*t (r')iy(«)(r, r')*„(r')*a(r) 



_ ^(m) _|_ ^^(m) _ 



= / dr 



^l{v)h^^\v)^^{v) + \j dr'*]„(r)*|„(r')W^("^)(r,r')*^(r')*n.(r) 



ly(am) ^ dr dr'4fl{r)^l{r')W^^'^\r,r')4f^{r')^a{r). 



(2) 
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The last term describes the conversion of atoms to molecules and vise versa and is given by 
[55, 57, 69, 76]: 



dr' 



r + r' 



^.{2a-^m)(r,r')*a(r)*„(r') + 



(3) 



The coordinates entering the field operators in (3) represent the annihilation (creation) of 

two atoms, one at position r the second at position r', and the creation (annihilation) of a 



molecule at the center-of-mass coordinate R = 
satisfy the usual commutation relations for bosons: 



The atomic, molecular field operators 



*„(r),*t(r') = *^(r),*J„(r') 



S{r-r') and *„(r),*„(r') 



= 0. Since the atoms and molecules are 

*„(r),#t(r') 



*„(r),*„(r') 

distinguishable, different particles, their mutual field operators commute, 

^a(r), ^'m(r') = 0. Finally, we note that the interaction terms appearing in the Hamilto- 
nian (1-3) are symmetric, i.e., W^^^r, r') = W^^'^r', r), . . . , #(™-2a)(^j.^ ^./^ ^ ^^(m^2a)(^j,/^ 
because the Hamiltonian is symmetric to the exchange of position of any two particles of 
the same kind. 

The Hamiltonian (1-3) commutes with the following particle- number operator 



N^Na + ^Nm^ Jdr [*t(r)#„(r) + 2*1„(r)*„(r)] , 



(4) 



refiecting a conservation law in presence of particle conversion. Accordingly, the Hilbert- 
space of the problem is a direct sum of Hilbert subspaces with different number of atoms 
and molecules: {N atoms; molecules}®{A'" — 2 atoms; 1 molecule}®{A^ — 4 atoms; 2 
molecules}® . . . ®{N — 2 [y] atoms; [y] molecules}, where [j] means the greatest integer 
not exceeding j. 

The purpose of this work is to treat the many-body Hamiltonian with atom-molecule con- 
version (1-3) multiconfigurationally. To this end, we expand the atomic ^'a(r) and molecular 
^^(r) field operators by two complete sets of time- dependent orbitals. 



*„(r) =^6fe(t)0fe(r,i). 



(5) 



The sets of atomic {0jfc(r, f)} and molecular {■^^/(r, i)} orbitals span the time- dependent 
Hilbert space in which the system is to be propagated. The advantages of time-dependent 
multiconfigurational expansions, see the Introduction, is the employment of time- dependent 
orbitals which change in time according to a time-dependent variational principle. This 
allows one to use in practical computations a smaller number of time- dependent orbitals 
than the number of time-independent orbitals that would have been required otherwise. A 
general multiconfigurational expansion, see section IV, employs M orbitals for the bosonic 
atoms and M' orbitals for the bosonic molecules. In particular, even if only one orbital 
is available for the bosonic atoms and another one for the bosonic molecules, the resulting 
theory goes beyond the Gross-Pitaevskii mean-field theory for this system [57], see subsequent 
section III. 

Finally, it is convenient to derive the relevant results first for the popular contact inter- 
action. 



iy(""^)(r, r') = XamS{r - r'), W^^''-^^\r, r') = W^^-^''\r, r') = X^J{r - r'). (6) 



Here and hereafter, the dependence of quantities on time is not shown explicitly whenever 
unambiguous. Below, we will work throughout sections III and IV A with the contact- 
interaction Hamiltonian (7) and handle the case of general interactions (1-3) thereafter, in 
section IV B. 



Thus, substituting Eqs. (5,6) into the generic Hamiltonian (1-3) we get 



^(2a^m) ^ J2 {<Pk \h^^^ I 0,) bib, + fYl MslMl) blblkK + 




(7) 
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III. THE SIMPLEST CASE OF ATOM-MOLECULE CONVERSION: CONVER- 
SION MEAN FIELD (FULLY- VARIATIONAL TWO-MODE APPROXIMATION) 



A. The multiconfigurational ansatz 

To introduce the nomenclature in the first stage of this work and, independently, as an 
interesting and relevant problem for itself, we consider the resulting theory when there is 
only one orbital available for the (bosonic) atoms and one orbital available for the (bosonic) 
molecules. The atomic orbital will be denoted by 0i(r, t) = 0a(r, t) and the molecular orbital 
by ipi{r,t) = ipm{r,t). The corresponding creation operators are denoted by b\{t) = bl{t) 
and c\{t) = cl^{t). The atomic and molecular creation, annihilation operators obey the 
bosonic commutation relations corresponding to the field operators. 

The problem we wish to solve may now be formulated. In the present section we would like 
to derive a multiconfigurational theory for atom-molecule conversion which is exact in the 
smallest Hilbert subspace possible for bosonic species, namely, the Hilbert space spanned 
by the single molecular orbital iprn{^,t) and single atomic orbital 0a(r, t). We term this 
specific case of the general theory: conversion mean field. More technically, this theory is a 
fuUy-variational extension of the literature two-mode approximation [56, 65] and, of course, 
of the Gross-Pitaevskii theory with conversion [57] . 

The multiconfigurational wavefunction takes on the following form: 

[7V/2] 

\^{t)) = Y,Cp{t)\N-2p,p;t), 

p=0 

1 , . N-2p 

1^ - = ^ , m)) {cl{t)y \vac) , (8) 

where \vac) is a common vacuum of no atoms and no molecules. The index p enumerates 
the number of bosonic molecules in the system. The corresponding number of atoms is 
— 2p. N is the maximal number of atoms in the system which is obtained when there 
are no molecules. Obviously, is an eigenfunction of the particle-number operator N, 

Eq. (4), with the eigenvalue N. The atomic and molecular number operators in the relevant 
Hilbert space boil down to Na = blpa and Nm = cl^Cm, respectively. The size of this Hilbert 
space is [y] -|- 1. 
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B. The functional action S of the time-dependent Schrodinger equation and its 
evaluation 



Solving the time-dependent Schrodinger equation with the Hamiltonian (7) and the mul- 
ticonfigurational ansatz (8) means finding the equations governing the time evolution of 
the atomic and molecular orbitals, 0a(r,t) and ■?/'m(r, i), and of the expansion coefficients 
{Cp{t)}. The derivation of these equations of motion for 0a(r, t), '0j„(r, t), and {Cp{t)} re- 
quires a time-dependent variational principle. We employ the Lagrangian formulation of the 
(Dirac-Prenkel) time-dependent variational principle [89, 90], also see Refs. [33, 53, 54], and 
write the functional action of the time-dependent Schrodinger equation which takes on the 
form: 



S[{Cp{t)},U^,t),7/;^{r,t)] 




dt 



-lia{t) [(0,(r,t)|0„(r,t)) - 1] - [(V^^(r,t)|V^„(r,t)) - 1] 

'[N/2\ 



-e{t) 



p=0 



(9) 



The time- dependent Lagrange multiplies f^a{t), fJ^m{t) and e(t) are introduced to ensure 
normalization of the atomic ^^(r, and molecular ip^{r,t) orbitals and of the expansion 
coefficients {Cp{t)} at all times. Hait) and Umit) also serve another role. They exactly "com- 
pensate" for those terms appearing within the Dirac-Prenkel formulation of the variational 



principle 



dt 



^(t) ) [20, 21], i.e., when the variation of ^(t) is performed 



before the expectation value {^{t) 



^{t)) is evaluated; see in this context 



[33, 90]. We shall see below and more elaborately in appendix Al that these Lagrange 
multipliers can be eliminated from the resulting equations of motion by making use of the 
normalization of the orbitals in combination with unitary transformations. 

The expectation value of the Hamiltonian in Eq. (9) can be expressed in two equivalent 
forms. The first form depends explicitly on the orbitals 0a(r,t), ?/'m(r, ^) and the second 
on the expansion coefficients {Cp{t)}. The two forms are needed to derive the respective 
equations of motion for the orbitals and expansion coefficients. 
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1. Orbital- explicit expression of S 



Utilizing the multiconfigurational expansion (8) and the individual terms of the many- 
body Hamiltonian (7), the first form of the expectation value of ^("^^) — i-^ reads: 



d 
di 



U"^ - i 



d_ 
di 



d_ 
di 



+ 

Ipm ) + 



+ - 



A. 



A. 



V2 



[N/2] 



p=0 



' dt 



(10) 



In Eq. (10) and hereafter we use the shorthand notation for expectation values of operators 

*(t)y etc. We 



with respect to *(t): (Na) = (^(t) ^{t)j, (blblc^j = (^^{t) blblcm 
can indeed see that expression (10) depends explicitly on the orbitals (f)a{r,t) and ipm{T',t) 
through integrals over one-body terms, two-body interaction terms, and the conversion term: 

(^(pa h}^"^ - i§i ^a)' {Mm \ Mm), etC. 

Making use of the multiconfigurational expansion (8), we can express the above expecta- 
tion values (^(f) |. . .| ^(t)) in a closed form. The expectation values of particle-conserving 
operators, hke the number operators Na and 7V^, with respect to ^{t) read: 

[N/2] [N/2] 

iV,) = J](iV-2p)|Cp(i)|% (7VJ = J]p|Cp(i)|% 



p=0 



p=0 



[N/2] 



Na{Na -1)) = ^{N- 2p){N - 2p - 1) \C,{t)f , 

p=0 

[N/2] [N/2] 

Nm{Nm - 1)) = E - 1) \^P(^)\' ' {NaNm) = ^ p(iV - 2p) \CM' " (H) 

p=0 p=0 

We see that the dependence of the expectation values (11) on the expansion coefficients 
is only through weighted sums Yl^^o^ of terms \Cp{t)f, i.e., that configurations of a 
different number p of molecules are not directly coupled. The expectation values of particle 
non-conserving operators, originating from the conversion of particles, are given by 

[N/2] 

[clkba) = Yl -2p+l){N-2p + 2) C;{t)Cp^i{t), 



bibicm 



"Irfiaba^^ i 



(12) 
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and seen to couple directly configurations with a different number of p and p — 1 molecules. 



2. Expansion- coefficient- explicit expression of S 



Utilizing the multiconfigurational expansion (8) and the many-body Hamiltonian (7) as 
a whole, we can express the expectation value of ^^(^a^m) _ jd_ functional action (9) 

as an explicit function of the expansion coefficients {Cp{t)}. One readily finds, 

d 



dt 



lN/2] 



lN/2] 

J2 (N-2p,p;t 



d 
di 



N -2p',p';t)Cpr -i 



dt 



(13) 



Eq. (13) contains yet another type of matrix elements, which are the representation of 
jj{2a±^m) _ -j-j^g subspace of configurations {\N — 2p,p]t)}. These matrix elements 

can be evaluated explicitly. We divide them into two types, recalling that the Hamiltonian 
is expressed as a sum of particle-conserving and particle non-conserving parts, — 
fyiam) ^ ^(2a^m)_ rpj^g diagonal, or particle -conserving matrix elements read: 



N-2p,p-t) = {N -2p){ct)a 



N-2p,p;t H^-^)-i^ 

at 

+^{N -2p){N -2p-l){ct>l\cl>l)+p 

+ ^P{P - 1) (0m I </'m> + KmP{N - 2p) {(j)a(t>m \ (t>a(t 



dt 



(pa 

UL 

p = 0,...,[7V/2]. (14) 



The off-diagonal, particle non-conserving matrix elements, originating from the conversion 
term W^(2«-"*) = wi^^^^) + wim^^a)^ ^.j^g following form: 



N - 2p, p; t 
A 



N-2{p-l),p-l;t 



y^p{N-2p + l){N-2p + 2){<l>m \<I>1), P = 1, . . . , [N/2], 



V2 

N-2{p-l),p-l;t 



yy{m^2a) 



N-2p,p;t 



p=l,...,[7V/2]. (15) 



N-2p,p;t^^ 
N-2{p-l),p-l;t)y , 

All other matrix elements of —i^m the subspace of configurations {\N — 2p,p;t)} 

vanish. With explicit expressions of the functional action (9) we can now proceed and derive 
the equations of motion of ^{t). 
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C. The equations of motion for '^{t) 



We perform the variation of the action functional (9,10) with respect to the orbitals 
and coefficients. Equating the variation of S [{Cp{t)}, (l)a{r, i), ■0m(r, t)] with respect to the 
orbitals to zero, eliminating the Lagrange multipliers /lait) and Umit) from the resulting 
equations (see appendix A 1 for details), and dividing the result by (^^a^ and (^Nrnj respec- 
tively, we obtained the following equations of motion for the orbitals: 



where the shorthand notation 0^ = -0^ = is used here and hereafter. The "inter- 
action strengths" are given by 



^Ar„(iV„-l)^ l^NaNra) {biblCm 



N„) (N„) (N, 



■Am(^) — Ajjj r , A^(j(/:) — ^am TZ. \ ) ■'^con(^) — ^con TZ. \ (1'^) 



and vary in time due to the conversion of atoms to molecules and vise versa. 

The quantities appearing on both the right- and left-hand sides of equations of motion 

(16) are projection operators and given by 

P(«) = 1 - , P("^) ^l-\'^m){i^m\. (18) 

When acting on one-body functions in the atomic and molecular spaces, P^"^ and P*^'") 
project these functions onto the subspaces orthogonal to the orbitals (f)a{r,t) and '0,„(r, i), 
respectively. The projection operators P*^"\ P*^"*) emerge when one ehminates the Lagrange 
multipliers fia{t), iJ^mif) from the equations of motion, see appendix Al. 

The appearance of the projection operator on both the left- and right-hand sides of 
Eq. (16) makes (16) a cumbersome set of two coupled integrodifferential non-linear equations. 
Can one simplify the matters? The answer is positive. 

To this end, we invoke the invariance properties of the many-particle wavefunction '^{t). 
Specifically, we can multiply the atomic 0(j(r, i) and molecular ip^{r,t) orbitals by time- 
dependent phase factors to give transformed orbitals 0„(r, i) and ip^{r,t), and from the 
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latter assemble transformed configurations \N — 2p,p;t) — > \N — 2p,p]t). Then, we can 
compensate for transforming the orbitals by the "reverse" transformation of the expansion 
coefficients {Cp{t)} {Cp{t)}. Overall, we write this invariance of the many-body wave- 
function as follows: 

W2] [N/2] _ 

\^{t)) = J2Cp{t)\N-2p,p-t)=J2Cp{t)\N-2p,p-t). (19) 

p=0 p=0 

Clearly, unitary transformations of the orbitals and the respective transformation of the 
expansion coefficients neither change the size of the Hilbert space nor couple configurations 
with a different number of molecules. To express these properties we use the same summation 
index p in both the middle part and right-hand side of Eq. (19). 

We can now make use of the invariance relation (19) to simplify Eq. (16). Specifically, 
there exists one unitary transformation that eliminates the projection operators acting on 
the time-derivatives (left-hand sides) in Eq. (16), without introducing any further constraint 
into the equations of motion; see appendix A 1 for more details. The equations of motion 
for the atomic and molecular orbitals thus finally read: 



Eq. (20) has the following property. Operating from the left with (0a | and (V'ml, respectively, 
we obtain the relations 

I k) = 0, ^V'm I V'm) = 0, (21) 

clearly ensuring that initially-normalized orbitals remain normalized for all times. We can 
see the meaning of the unitary transformation carrying Eq. (16) to Eq. (20). This unitary 
transformation takes normalized time-dependent orbitals, (0a | 0a) = 1 and ('0^ | ■0^) = 
1, which therefore satisfy the general relations ^^''''"['^"^ = ^0a I 0a^ + ^0a I 0a^ — and 
d{n>,-^ti>m) _ I _|_ I -0^^ = 0, and transforms them to time-dependent orbitals 
satisfying the specific differential condition (21). 

Before we move to the corresponding working equations for the expansion coefficients 
{Cp(t)}, it is instructive to enquire whether we could further simplify the equations of 
motion (20), by eliminating the projection operators P^"\ p("*) also from the right-hand 
sides. The answer is in general negative. If we could eliminate the projection operators 
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remaining on the right-hand sides, it means that conditions (21) are not satisfied any more. 
What would then guarantee that the atomic and molecular orbitals remain normalized at 
all times? It turns out that the condition for that is: Im {Aeon (c^^a^a) (^m I ^a^a)} = for 
all times. In turn, even if this condition is satisfied at t = 0, it is not in general guaranteed 
that it remains so for all times. Thus, the presence of particle conversion does not allow 
one to eliminate the projection operators P*^"\ p("*) also from the right-hand sides of the 
equations of motion (20). Alternatively speaking, in the absence of particle conversion, it is 
possible to eliminate the projection operators completely from (20), see in this context [91]. 

To derive the equations of motion of {Cp(t)}, we equate the variation of the action func- 
tional (9,13) with respect to the expansion coefficients to zero and eliminate the Lagrange 
multiplier e{t) (see for details appendix A 1). The following equations of motion are obtained: 



?i(2a^™)(t)C(t)=Z^^, 

n^Zr'^\t) ^(n- 2p,p; t ^(^«^-) - 4 ^ - ^P'^P'-^ ' (22) 

where the vector C{t) collects the expansion coefficients {Cp{t)}. Eq. (22) is a set of cou- 
pled first-order differential equations with time-dependent coefficients, and preserves the 
norm of an initially- normalized vector of coefficients C(0). The time-dependent coefficients 
Ti!pp7^^\t)^ being the matrix representation of H _ in the subspace of configu- 

rations {|A^ — 2p,p;t)} and hence depending on the atomic (j)a{T^,t) and molecular '^„i(r, t) 
orbitals, are prescribed in the previous subsection lllB. 

Next, we make use of the invariance of the multiconfigurational wavefunction to unitary 
transformations (19). Explicitly, the unitary transformation responsible for transforming 
Eq. (16) for the orbitals to Eq. (20), transforms equations of motion (22) for the expansion 
coefficients to the final form: 

jj(2a^m)(^)C(i)=i^^, 

HiT"^^ {t) = (N- 2p, p; t (2'^--) N - 2p', p'; t) . (23) 



Equivalently, Eq. (23) can be obtained from Eq. (22) by substituting into the latter the 
differential condition (21). 

The coupled sets of equations of motion for the atomic 4>{r,t) and molecular '0^(r, i) 
orbitals and expansion coefficients {Cp{t)}, Eqs. (16) and (22) or, respectively, Eqs. (20) and 
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(23) constitute the conversion mean field theory (fully- variational two- mode approximation) 
for the interacting atomic-molecular system with conversion. 



D. The stationary self-consistent coherent mean field (time-independent fuUy- 
variational two-mode approximation) 

The theory presented above is a time-dependent many-body theory. It is certainly inter- 
esting to enquire what are the corresponding stationary many-body states of the atomic- 
molecular Hamiltonian (1-3)? In other words, what are the self-consistent solutions that 



minimize (extremize) the expectation value yi/ 

.[N/2] 



^) for a given time-independent 



multiconfigurational ansatz |^) = ^pJo Cp\N — 2p,p) assembled from time-independent 
atomic 0a (r) and molecular orbitals ■?/'m(r)? To get this stationary self-consistent coher- 
ent mean field (time- independent fully- variational two-mode approximation), we resort to 
imaginary time-propagation and set t — > —it in the corresponding equations of motion. 

Setting t — > —it in (the left-hand side of) either Eq. (16) or (20), the left-hand side 
decays to zero in time and the equation becomes time- independent. Then, by multiplying 
the result, respectively, by (^^aj and (^Nm^, and translating back the projection operators 
P(") and P('') to the corresponding Lagrange multipliers fia and fit, (see appendix Al), we 
obtain the multiconfigurational self- consistent (time-independent) equations for the atomic 
and molecular orbitals: 



Na) h^"-^ + A„ {Na{Na - 1)) \(t>a? + Km {NaNm) \lpm\'\ \<Pa) + 

+V2Xco„ (blblcm^ 0* \lprn) = l^a |0a) , 
nJ) + (iV^(iV^ - 1)) \lPm\' + Km (NaNm) \<Pa\'] H'm) + 

+ ^ (^clbaba^ (j)a\(l>a) ^ f^m\lpm) ■ (24) 

Similarly, restoring the Lagrange multiplier e{t) into either Eq. (22) or (23), see in this 
respect appendix A 1, and setting t — > —it therein, we obtain the stationary (self-consistent) 
eigenvalue equation 



jj(2a^)C = sC, 



N-2p',p' 



(25) 
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for the expansion coefficients. We see that the (redundant) time-dependent Lagrange multi- 



pher e{t) of the time-dependent theory has emerged as the eigenenergy e = (^'^ 
of the stationary theory. 

The theory "distilled" into Eqs. (24,25) is a fully-variational stationary theory for the 
interacting atomic-molecular system in presence of conversion, where a single orbital is 
allowed for the atoms and a single orbital to the molecules. It is a system of coupled 
eigenvalue-like equations for the orbitals and eigenvalue equation for the coefficients, thought 
non-linear and integrodifferential ones. 



IV. THE GENERAL MULTICONFIGURATIONAL THEORY WITH ATOM- 
MOLECULE CONVERSION 

In this section we develop a general many-body theory for atom-molecule conversion, by 
allowing the atoms and molecules to occupy more orbitals. Section IV A builds the theory 
for the popular contact interaction, whereas the case of generic non-contact interactions is 
presented in section IV B. 



A. Formulation for contact interactions 

1. The multiconfigurational ansatz for the wavefunction 

The multiconfigurational expansion mixes atomic-molecular states with different numbers 
of particles which are eigenfunctions of the particle- number operator TV |^(i)) = |^(i)): 

lN/2] 

\f}F,wF]t) = X 

X {h\{t)f' ■ ■ ■ (6L(t))"" {Ait))"'' . . . (f^,{t)Y^' \vac) . (26) 

We collect the individual occupations in the vectors fV = (nf , . . . , n^), mP — (mi, . . . , m^,). 
The number of bosonic atoms = n{ + . . . + n\j — N — 2p and molecules |m^| = 
nil + ■ ■ ■ + TTi^, — p oi each configuration \nP, mP; t) satisfies the particle-conservation law 
|fi^| -I- 2|77?.*'| = N. Observe that the number of molecules p serves as an index to the 
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occupation numbers fP and rh^. This simply reflects the fact that, for a given number 
oi N — 2p atoms and p molecules, the possible occupation numbers which the configu- 
rations can assume depend on p itself. The index p together with the occupation num- 
bers n^, mP make a unique representation of each configuration. The atomic and molec- 
ular number operators in the corresponding Hilbert space boil down to Na = Yl^=i ^I^fc 
and Nrn — 'Y^w=i^k'^k'i respectively. The size of the resulting Hilbert space is given by 

I N -2p + M -l\ fp + M' -l\ 
z2 'n h i-^-' by the sum of products of the sizes of the re- 

^ y M-1 j y M'-i j 

spective Hilbert subspaces for N — 2p bosonic atoms with M orbitals and p bosonic molecules 
with M' orbitals. 



2. Reduced density matrices for systems with particle conversion 

As part of the variational derivation we will need the expectation value of — 
with respect to ^'(t). To this end, it will be proved valuable to define and employ the 
reduced density matrices of ^(i). We remind that Lowdin has introduced the concept of 
reduced density matrices for systems of a fixed number of particles (identical fermions) [44] . 
Nevertheless and although ^(i) is not comprised of a fixed number of atoms or a fixed 
number of molecules, it is possible to define the reduced density matrices of a mixture of 
atoms and molecules with conversion. 

Having at hand the normalized many-body wavefunction ^(t), the reduced one-body 
density matrices of the atoms and molecules are defined by: 

M 

p(«)(ri|r2;i) = (*I(r2)*.(rO) = J2 P^S{t)<t>lM<i>,{v,,tl 

k,q=l 
M' 

p(-)(ri|r2;t) = {^l{r,)^m{r,)) = J] P^^^^^mi'Mi^^'M, (27) 

k',q'=l 

where the matrix elements p^^qi^) ~ (Hfiq^ ptjq'i^) ~ (^^I'^g'^ prescribed in ap- 
pendix B. We collect these matrix elements as p*-"-'(i) = \ Pk^{t) \ and p^"^\t) — I p^^qi{t) \. 
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Similarly, the reduced two-body density matrices of the atoms and molecules are defined by: 
p(«)(ri,r2|r3,r4;i) = (*l(r3)*l(r4)*„(r2)*a(ri)' 



M 



= E P?iW'^fe(r3,t)C(r4,t)0Kr2,t)0,(ri,t), 

k,s,l,q=l 

p(™)(ri,r2|r3,r4;t) = (#J„(r3)*]„(r4)*^(r2)*^(ri)) = 

M' 



k',s',l',q'=l 



P 



[am) , 



(ri,r2|r3,r4;t) = ( *i(r3)*a(ri)*L(r4)*n^(r2) 



M M' 



= E E pK'(^)'^^(r3,^)0.(ri,t)^:'(r4,t)^,'(r2,t), (28) 

fc,?=l fe',q'=l 

where the matrix elements = (blblkbg^, pt^^'vA^) = (cI'CyQ/c^/^, and pt^^q,{t) = 

^6^6gc[,Cg/^ are prescribed in appendix B. Because the reduced density matrices (27) and 
(28) directly only couple configurations with the same number of atoms and molecules, 
we will refer to them as particle- conserving reduced density matrices. In this context, 
p^""^^(ri, r2|r3, r4; t) is the lowest-order inter-species particle-conserving reduced density ma- 
trix. 

Prom the above discussion it is anticipated that, due to the conversion term (3) in the 
Hamiltonian, another kind of reduced density matrices appear in the theory. Specifically, 
we define the particle non- conserving reduced density matrices as follows: 

M' M 

p(2--)(r„r2|r3;t) = {^l{v,)i^ a{v2)^ a{r,)) = E E pg"H^)V':'(r3, ^)0.(r2, t)0,(ri, t), 

k'=l k,q=l 
M' M 

^'^)(r3|r2,ri;t) = (*i(ri)*i(r2)*^(r3)) = E E 4r/7"^^(^)'/';(ri, ^)0:(r2, t)^.'(r3, t), 

k'=l k,q=l 

p(— 2«)(r3|r2,r,;t) = {p(^"\r„r,\r,;t)y , p%7''\t) = {pg^^^t)}*. (29) 
The matrix elements pf,^^ = ((^kihkbq) are given in appendix B. 
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3. The functional action S and its evaluation 



We start from the functional action of the time-dependent Schrodinger equation which 
in the general multiconfigurational case takes on the form: 



S[{Cnr>Mt)},{Mr,t)},{'ilJk'(r,t)}]^ J dtl (^^{t) 



dt 



M 



M' 



k,j=l k',j'=l 
'[N/2] 



-S{t) 



p=0 nP,mP 



(30) 



The time-dependent Lagrange multiplies {fJ'kjit)}: {l^k'^it}} ^{t) introduced to 
ensure orthonormalization of the atomic {0fc(r,t)} and molecular {tjjkilrjt)} orbital sets and 
normalization of the expansion coefficients {C^Pmp(^)}- 

To derive the equations of motion for the atomic-molecular multiconfigurational wave- 
function (26), the expectation value of "i) _ -^th respect to ^{t) is needed, where 
fj{2a^m) -g gjygji gq (jy 'pjjg expectation value of _ jg expressed by two 

equivalent forms, as done in section IIIB. The first form, where the dependence of Eq. (30) 
on the atomic and molecular orbitals is explicit, reads: 



dt 



\ ^ I 

' fc,o=l \ 



M 



•',1 
M' 



- i — 



(rn) 



+ % Yl Pti'l'q' {^k'(l>s'\((>q'((>l') + ^am ('^fc^fc' I <^ A' ) + 



k,s,l,q=l 
M' 



k',q'=l 



dt 



d 



+ 



^a' ) + 



(31) 



M M' 



k',s',l',q'=l 

M M' 



, m ivi 

yPk'kq Wk' I Mq) + P\kk 
^ k,q=lk'=l 



k,q=l k',q'=l 

,^2a) 



{4>q4>k I V'fe' 



[N/2] 



-EEC 

p=0 nP,faP 



dCfiPjfiP 



We see in (31) the appearance of the particle-conserving and particle non-conserving reduced 
density matrices introduced in the previous subsection IV A 2. Eq. (31) is to be used to derive 
the equations of motion of {(f)k{v,t)} and {■0jt'(r, t)}. 
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The second form of the expectation value of _ in the functional action (30), 



dt 



[JV/2] 

p=Q fiP,mP 



X 



X 



[N/2] 

p'=0 n'p\m'p' 



d_ 
di 



->/p' -ip' \ ri -^^^ 



nPfnP 



dt 



(32) 



displays its explicit dependence on the expansion coefficients, and therefore will be employed 
to derive the equations of motion of {Cf;p^p(t)}. Finally, it is deductive to compare the 
structure of Eqs. (10,13) in the conversion mean field (fuUy-variational two-mode) problem 
to that of Eqs. (31,32) of the general problem. 



4- The equations of motion for ^(i) 



Collecting the above ingredients, wc arc ready to perform the variation of the functional 
action 5* [{Cjjp^p(t)}, {^^(r, t)} , {■^^/(r, t)}] and arrive at the equations of motion of ^(i). 
Equating the variation of the action functional (30,31) with respect to the orbitals to zero 
and eliminating the Lagrange multipliers {A*fe"^(^)}, {l^k^j'it)} (see appendix A 2), we obtain 
the following result, j = 1, . . . , M, / = 1, . . . , M': 



M 



M 



M' 



k=l 



q=l 



k'=l 



M' / W , M 

k'=\ \o'=l ^ k,q=l 



, (33) 



where terms with products of reduced two-body density matrices times orbital pairs are 
collected together and denoted for brevity as 



M 



M' 



{p2{<P',i^')Y:^^Kj2pklqim)+Krn E P^^qq'i^l'M 



!,l = l 

M' 



k',q'=l 



M 



s',l'=l k,q=l 



(34) 
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and 

M M' 

P(«) = 1 - ^ , pM = 1 - 5^ [ilJu' ){M , (35) 

u=l u'=l 

are projection operators. When acting on one-body functions in the atomic and molecular 
spaces, P(") and P*^"^) project them onto the subspaces orthogonal to those spanned by the 
orbitals {0A;(r, t)} and {ipk'ir, t)}, respectively. These projection operators emerge when one 
eliminates the Lagrange multipliers {fJ'^^jit)} and {Ai[™^(t)} from the equations of motion, 
see appendix A 2. To remind, we use the shorthand notation = ipji = in the 
equations of motion for the orbitals. 

The appearance of the projection operators on both the left- and right-hand sides of 
Eq. (33) makes the system (33) a cumbersome system of integrodifferential non-linear equa- 
tions. This situation can be simplified by generalizing the treatment of section III C, namely, 
exploiting the invariance properties of the many-particle wavefunction ^(t). Specifically, we 
perform independent unitary transformations on the atomic {0fc(r,t)} — |0^(r,t)} and 
molecular {'^^/(r, t)} —>■ {^ipi^,{r,t)} orbitals, which results in transformed configurations 
\fP','mF;t) — > \fiP,mP;t). Then, we can compensate for the transformations of the orbital 
sets by the "reverse" transformation of the expansion coefficients {CjiPrnp{t)} — > {C'fjp^(i)}. 
We represent this invariance by the following equality: 

= E E Cnr'Mt) |n^mf;t) = E E CH.rnp{t)\nP,mP;t). (36) 

The transformations of the orbital sets and expansion coefficients do not change the size of 
the Hilbert space, or couple systems with different numbers of atoms, molecules. We remark 
that transformations which inter-mix atomic and molecular orbitals are not required for our 
needs. To represent these properties, the same occupation numbers n^, and summation 
index p of the number of molecules are used for both multiconfigurational expansions of ^(t) 
in Eq. (36). 

We can now make use of the invariance (36) to simplify the equations of motion (33), 
without introducing further constraints into the equations of motion. We utilize a specific 
unitary transformation of the many-particle wavefunction that eliminates the projection 
operators acting on the time-derivatives (left-hand sides) in Eq. (33); see appendix A 2 for 
more details. The final result for the equations of motion for the atomic {0jfc(r, i)} and 
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molecular {■0fe'(i", t)} orbitals thus takes on the form, j = 1, . . . , M, / = 1, . . . , M': 



i Uj 



p(a) 



M 



M 



M' 



+E {/^^"H^)};' X y: {P2(0^ ^^)}l;^ 10.) + V2x,on Yl i^.') 



k=l 



p(m) 



M' 



k'=l \q'=l ^ k,q=l 

with the projection operators P*^"^ and P^'") appearing now on the right-hand sides only. 

Now, taking the respective scalar products of (37) with {(0jt|} and {(V'jk'l}) '^^ obtain the 
following differential conditions: 



k'=l 



M 



(37) 



(f>k I 0? 



0, k,q^l, 



,M, 



V'fe' I i^q 



0, k',q' = l, 



M'. 



(38) 



It is instructive to mention that these differential condition have been introduced originally 
by the MCTDH developers [7, 8], and used thereafter in particle-conserving multiconfigura- 
tional theories for identical particles and mixtures [29-33, 53, 54]. 

The differential conditions (38) ensure that initially-orthonormalized orbital sets 
{0fe(r, t)}, {■0fc/(r, i)} remain orthonormalized at all times. The meaning of the uni- 
tary transformation carrying equations of motion (33) to (37) can now be seen. This 
unitary transformation takes orthonormal time-dependent orbitals, {(pk \ <f>q) = Skq and 
{'ipk' I V'g') = ^k'q' [^kq, ^k'q' is the Dirac delta-function], which therefore satisfy the general 



relations ^^^^ 



4>k I 0« > + ( 0fe I 0, 



and 



= {^k' I ^q) + {^k' I ^q' 



0, 



and transforms them to time-dependent orbitals satisfying the specific differential conditions 
(38). 

Moving to the equations of motion for the coefficients {Cfjp^p(i)}, we equate the variation 
of the functional action (30,32) with respect to the expansion coefficients to zero. Eliminating 

the Lagrange multiplier e{t) by a respective phase transformation of the coefficients, we arrive 
at the form: 



nPmP ,n'p' m'p' ^ ' 



nF, wF\ t 



at 



(39) 
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Eq. (39) has exactly the same form as Eq. (22) in the specific theory of section IIIC, and 
constitutes a set of coupled first-order differential equations with time-dependent coefficients 
that preserve the norm of an initially-normalized vector of coefficients C(0). The matrix 
elements of ^(2«— H _ {A with respect to two general confierurations, Tii^"^™^, , ^, At), are 
prescribed in appendix B. 

Now, to arrive at the final form of the equations of motion for the coefficients, we make 
use of the invariance of ^(i) to unitary transformations (36). Specifically, the unitary 
transformation carrying Eq. (33) for the orbitals to Eq. (37), casts Eq. (39) for the expansion 
coefficients into the final result: 

H(2a^m)(^)C(i)=i^2W, 

^W-'^v it) = {f^,rn';t\ if (^"--) I n'^'m'^' ; t) . (40) 

Eq. (40) has exactly the same form as Eq. (23) in the specific case of section III C. Of course, 
there are much more expansion coefficients in the general case. Finally and equivalently, we 
note that the result (40) can be obtained from Eq. (39) when the differential condition (38) 
is substituted into the latter. 

Let us pause for a moment and summarize. We have started from the functional action 
(30) and arrived at the equations of motion for ^{t). Eq. (33) for the orbitals {0fc(r, i)}, 
{ipk'{Y,t)} and Eq. (39) for the expansion coefficients {CrtPmpit)}, or, respectively, Eqs. (37) 
and (40) constitute the multiconfigurational time-dependent Hartrcc theory for systems with 
particle conversion; here specifically the theory for bosonic atoms and bosonic molecules with 
conversion. 



5. The stationary self-consistent general multiconfigurational theory with conversion 

The theory presented above is a time-dependent many-body theory and, as done in the 
previous section III, it is relevant to put forward the corresponding stationary general theory. 
Consider the multiconfigurational expansion |\E') = Y2^^o^ '^hp mp ^nPmP \TP',rnJ'), where the 
expansion coefficients {C^jp^p} and orbitals {0fc(r)}, {■^^/(r)} assembling the configurations 
{|ri^, TvP)} are time- independent quantities. What are then the self-consistent solutions that 
minimize (extremize) the expectation value ( ^ 
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The working equations of the stationary theory can be obtained by resorting to imagi- 
nary time-propagation and setting t — > —it in either equations of motion (33) and (39) or 
equations of motion (37) and (40), after the time-dependent Lagrange multipher s{t) has 
been reinstated. Then, by translating the projection operators P'^"\ pM to the respective 
Lagrange multiphers {a*^"'*}, {l^k^j'}, the resulting working equations take on the form: 



M 

E 

q=l 
M' 

E 

q' = l 



p^>^"^ + {P2(0^v^^)}l;^ 

k = l,...,M, 

p^/^^"^^+{P2(^^<^^)}^?j 



M' 



M 



^con 

72 



fe'=l 



M 



M' 



i'=i 



fc,o=l 



(41) 



for the orbitals, and 



(2a^m) 



^ (2a^m) 



(42) 



nPmP M'p' m'p' \ 

for the expansion coefficients. As seen in section III, the time-dependent Lagrange multipher 
e{t) emerges in the time-independent theory as the eigenenergy of the coupled atom-molecule 



system with e — 



The stationary equations for the orbitals (41) can be further simplified. The stationary 
wavcfunction \[', as its time-dependent counterpart, is invariant to independent unitary 
transformations of the orbital sets {0i;(r)}, {'^^'(r)} and the "inverse" transformation of 
the expansion coefficients {CfiPij^}. We can use the unitary matrices which diagonalize the 
matrices of Lagrange multiphers {jJ^^^]}-, {A*fey}- We note that the matrices of Lagrange 
multipliers are Hermitian matrices for stationary states. As a result of this transformation, 
we obtain a set of coupled equations for the orbitals that look just as Eq. (41), except for 
the right-hand sides being diagonal: 



M 

E 

q=l 

M' 

E 



A; = 1, . . . , M, 

pS/^(-) + {p,(v^^0^)}l?j 

k' ^l,...,M'. 



M' 



k'=l 



M 



V2 



k,g=l 



(43) 
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The form of the equation for the expansion coefficients (42) does not change. Thus, the 
final result for the stationary theory, Eqs. (42,43), is a coupled system of integrodifferential, 
non-linear equations constituting eigenvalue-like equations for the orbitals and eigenvalue 
equation for the expansion coefficients; compare to section HID. Generally, the transfor- 
mation of the matrices of Lagrange multipliers to diagonal form would make the orbitals 
delocalized. Hence, in problems where working with localized orbitals is of advantage or rel- 
evance, for instance in lattices, it is the form (41) with the in-general non-diagonal Lagrange 
multipliers which is to be preferred. 



B. Formulation for general interactions 



The last stage of the theory is to return to the case of generic non-contact interactions 
in the Hamiltonian (1-3) and have the respective theory derived. The derivation of the 
equations of motion follows essentially the same steps taken in the previous subsection IV A 
and there is obviously no need to repeat it. The only extra care needed is when minimizing 
the conversion term in the functional action with respect to the molecular orbitals, where 
exchange of variables is used. This point and related derivations are discussed in appendix 
A 2. Below, wc report the final results of the time-dependent as well as the self-consistent 
time-independent theories. 

The final form of the time-dependent equations of motion of the orbitals reads, j — 
1,...,M,/ = 1,...,M': 



m4- ) = p^"^ 



M 



M 



M' 



k=l 



k'=l 



M' / M' M 

k'=l \g'=l ^ k,q=l 



(44) 



(2a^m) TT r(2a^m) 

k'=l \q'=l ^ k,q= 

where terms with products of reduced two-body density matrices times one-body potentials 
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(see below) are collected together and denoted for brevity as 

M M' 

s,l=l k',q'=l 
M' M 

- E + E pS'<"^- (45) 

s',l'=\ k,q=l 

Comparing Eq. (37) of the previous subsection IV A to Eq. (44), we see that in the latter 
more general, time-dependent local potentials appear which are given explicitly by: 

W^f{r,t)= J cl>:{r',t)W^^\r,r')^i{v',t)dv', 
^5Hr,t)= J rA^',t)W^"^\T,r')ij,{T',t)dr', 

Wtr\^,t)- J cf>l{r',t)W^'''^\v,r')cf>,{r',t)dr', 

<~^(r, t)=J dr'W^^^--^ (^r + ^, r - 0, + 0, - ^, , 

W^ir'"^(r,t) = / rfr>;(r',t)W^(— ^«)(r,r')^.^(^,t) • (46) 

These potentials derive from the interaction terms and conversion term in the Hamilto- 
nian (1-3) and, in the specific case of contact particle-particle interactions, boil down to 
products-of-orbitals, see for comparison Eq. (37). The form of the equations of motion for 
the corresponding expansion coefficients, 

does not change for general interactions. Of course, the matrix elements i?i^°i~'"!L do 
depend on the specific form of the particle-particle interactions. 

Finally, the self-consistent, time-independent general theory is obtained from the time- 
dependent one by taking t — > —it. The stationary self-consistent equations for the orbitals 
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read: 

M 

E 

<?=i 



M' 



M 



{ra—^2a) jtr (m^2a) 



k'=l 



,M, 



M' 

E 



k=l, 
k' ^1,...,M', 



fe,o=l 



/=1 



(48) 



where the arrows indicate the Lagrange muhiphers in their diagonal form, as done in 
Eq. (43). Finally, the self- consistent-eigenvalue form of the equation for the expansion 
coefficients. 



j^{2a^m)Q = eC 



nPmP ,n'p' m'p' 



(49) 



remains unchanged for general interactions. 

Eqs. (44) and (47) constitute a multiconfigurational time- dependent theory for systems 
of bosonic atoms and molecules with conversion (particle conversion in the generic case). 
Furthermore, Eqs. (48) and (49) constitute a multiconfigurational self-consistent time- 
independent theory for systems of bosonic atoms and molecules with conversion (particle 
conversion in the generic case). Both theories extend the scope of the successful multi- 
configurational time-dependent Hartree method and its versions specified for systems of 
identical particles and mixtures to new physical systems and problems. 



V. SUMMARY AND CONCLUDING REMARKS 



In this work we have derived a many-body propagation theory for systems with particle 
conversion. The theory is intended for systems with a finite number of interacting particles, 
typically in a trap potential. The theory has been exemplified and working equations have 
been explicitly derived for systems of interacting structureless bosonic atoms and bosonic 
molecules undergoing the conversion 'reaction' 2a ^ m. In doing so, we have also extended 
the scope of the successful multiconfigurational time-dependent Hartree (MCTDH) method 
and its versions specified for systems of identical particles and mixtures to new physical 
systems and problems. We note that the MCTDH method is considered at present the most 
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efficient wave-packet propagation approadi for in-general distinguisfiable coupled degrees- 
of-freedom, with no particle conversion, of course. The general multiconfigurational theory 
with particle conversion shall be referred to as MCTDR- conversion, whereas the explicit 
scenario derived throughout this work by MCTDH-[2a^m]. 

To treat systems with particle conversion, one has to work in second quantization formal- 
ism, where the Hamiltonian with particle-conversion terms can be represented. The next 
step is to define the configurations. In presence of particle conversion configurations with 
different numbers of atoms and different numbers of molecules are coupled. For instance, 
consider the particular case of the conversion 'reaction' 2a ^ m. In this case, the subspace 
of coupled configurations can be easily obtained by starting from the configurations made of 
N atoms only, and operating repeatedly with the conversion operators in the Hamiltonian 
until configurations made of the maximal number [y] of molecules are reached. 

In the multiconfigurational theory for the 'reaction' 2a ^ m, there are M time- dependent 
orbitals {(f)k{r,t)} available for the atoms and M' time- dependent orbitals {■0fe'(r, i)} for 
the molecules. The multiconfigurational ansatz for the many-particle wavefunction ^(i) 
is taken as linear combination with time- dependent coefficients {C^Pmp(i)} of all possible 
configurations {|n^, m^; t)} assembled from p molecules and N — 2p atoms - distributed over 
the M and M' respective orbitals - and coupled by the conversion term in the Hamiltonian. 

The evolution of ^(t) is then determined by the Dirac-Prenkel time-dependent variational 
principle. Utilizing the Lagrangian formulation of the Dirac-Prenkel variational principle, 
one arrives at two sets of coupled equations of motion: The first set is for the orbitals 
{0fc(r,t)} and {^/'^/(r, t)}, and the second for the expansion coefficients {C^Pmp(^)}- The 
first set is comprised of first-order differential equations in time and non-linear intcgrodiffer- 
ential (for non-contact interactions) in position space. The second set consists of first-order 
differential equations with coefficients forming a time-dependent Hermitian matrix. Thus, 
equations of motion (33,39), or Eqs. (37,40), constitute the time-dependent multiconfigura- 
tional theory for bosonic atoms and molecules with conversion - MCTDH-[2a^m]. 

The structure of the equations of motion for systems with particle conversion reminds 
of the structure of the equations of motion in multiconfigurational time-dependent theories 
for systems of identical particles and mixtures [32, 33, 53, 54]: (i) There are projection 
operators on the right-hand-sides of the equations of motion for the orbitals, ensuring that 
the respective orbitals remain normalized and orthogonal to one another for all times; (ii) 
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The equations for the expansion coefficients are first-order differential equations with time- 
dependent coefficients; and (iii) The equations of motion for the orbitals are formulated 
in terms of reduced density matrices. This resemblance would allow one to transfer the 
effective numerical techniques that have been developed in the past almost-twenty years for 
multiconfigurational time-dependent many-body systems without particle conversion [7, 8, 
12, 29-33, 53, 54] to the present theory for systems with particle conversion. 

Particular attention has been paid to the reduced density matrices appearing in the 
theory. As the multiconfigurational expansion involves configurations with different numbers 
of atoms and molecules, two types of reduced density matrices are defined. There are 
particle- conserving reduced density matrices, p*^")(ri|r2; t), p^"^\ri\r2;t), p*^"''(ri, r2|r3, t), 
p^"^\ri, r2|r3, r^; t) and p*^"™')(ri, r2|r3, t), which directly do not couple configurations with 
different numbers of atoms and molecules. Despite this property, the particle-conserving 
reduced density matrices, such as p^°'\ri\r2;t) and p(")(ri, r2|r3, i), are no^ the standard 
density matrices introduced by Lowdin [44] for many-particle systems without conversion. 
The second type of reduced density matrices that appear in the theory are particle non- 
conserving reduced density matrices, p^°'^'^\ri, r2\r^; t) and p^"*^"^(r3|r2, fi; t), and originate 
from the conversion term in the Hamiltonian. They obviously have no analogs in systems 
without conversion. Here, it is of interest by itself to study properties of particle-conserving 
reduced density matrices and certainly of particle non-conserving reduced density matrices 
in systems with particle conversion. 

The time-dependent multiconfigurational theory MCTDH-[2a^m] readily admits the cor- 
responding stationary theory. By resorting to imaginary time propagation, the equations of 
motion of the time-dependent theory boil down to the fully-self-consistent time-independent 
multiconfigurational theory, Eqs. (48,49), for stationary states of the system 2a ^ m, in pres- 
ence of all particle-particle interactions, of course. With this result, available self-consistent 
multiconfigurational theories for systems without particle conversion, noticeably for fermions 
[92, 93], distinguishable degrees-of- freedom [13], bosons [52], and mixtures [54], are taken a 
step further, to systems with particle conversion. 

A specific case of interest for systems of bosonic atoms and molecules with conversion 
is the case of M = 1 atomic and M' = 1 molecular orbitals, which is presented in section 
III before the general MCTDH-[2a^m] theory is developed. For M — 1 and M' — 1, 
the corresponding multiconfigurational theory is the fuUy-variational theory that results 
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when the shape of the atomic 0a(r, t) and molecular ipmi'i'^, t) orbitals and of each and every 
expansion coefficient Cp{t) are optimized according to the variational principle. Being fuUy- 
variational with respect to the shape of the orbitals '?/'m(r, t) and with respect to 

the expansion coefficients {Cp{t)}, the theory generalizes the literature Gross-Pitaevskii 
equation [57] and two-mode approximation [56, 65] for bosonic atoms and molecules with 
conversion. We term this specific case of the general theory conversion mean field, as there 
is only one orbital available for the bosonic atoms and one for the molecules - the minimal 
number possible for bosonic species. 

At the other end, in the limit where the number M of atomic orbitals and M' of molecular 
orbitals goes to infinity, the MCTDH-[2a^m] theory becomes an exact representation of 
the timc-dcpcndcnt many-particle Schrodinger equation with the particle-conversion Hamil- 
tonian ^(2a^"i) practice, one obviously has to limit M and M'. Here, the employment 
of time- dependent orbitals, which has been very successful for the MCTDH approach and 
its versions specified for identical particles and mixtures, is of great help and advantage. Of 
course, even with time-dependent orbitals the size of the Hilbert space grows rapidly with 
the size of the system and the number of orbitals M and M' employed. Consequently, with 
increasing system size and as the number of orbitals which one has to employ becomes larger, 
e.g., for stronger interactions, it is instructive to devise truncation schemes beyond the us- 
age of time- dependent multiconfigurational expansions over complete Hilbert subspaces. We 
mention two such truncation strategies: (i) to truncate time-dependent multiconfigurational 
expansions to include parts of Hilbert subspaces, i.e., to include not all available configura- 
tions for a given system size and number of orbitals M,M'; and (ii) to concentrate on the 
reduced density matrices, write equations of motion for them directly, and thereafter trun- 
cate the resulting hierarchy of equations of motion for higher-order reduced density matrices 
at some given order. The development of these truncation schemes for time- dependent mul- 
ticonfigurational expansions in systems with particle conversion extends beyond the scope 
of the present work. 

Finally, the explicit equations of motion presented in this work are for the specific 'reac- 
tion' 2a ^ m where the atoms and molecules are structureless bosons - the MCTDH- [2a^m] 
theory. Several other systems come to mind: (i) Other 'reactions' with bosonic atoms of 
the same kind, e.g., 3a ^ m; (ii) 'Reactions' with bosonic atoms of a different kind, e.g., 
a -\- a' ^ m. In this case and for general particle-particle interactions, the center-of-mass 
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coordinate is, of course, R = ^mXm7 ' '^^^^^ "^a' are the masses of the respective 

species; (iii) 'Reactions' including fermionic atoms, e.g., o/ + a?, ^ m/ and aj + ^ rrib 
where the subscript b, f stands for bosonic, fermionic species. In the latter case, a unified 
form of the respective equations of motion and those of the present work is anticipated; and 
(iv) A whole zoo of 'reactions' for particles with spin, internal-structure. The extension of 
MCTDH-[2a^m] for the above concrete examples as well as for other systems with particle 
conversion can be done by following the theory and derivation steps of the present work. 
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APPENDIX A: FURTHER DETAILS OF THE DERIVATION OF THE EQUA- 
TIONS OF MOTION 



1. The specific case of the conversion mean field (fully- variational two- mode ap- 
proximation) 

When the variations of the functional action (9) with respect to the orbitals and expansion 
coefficients are put to zero, the following equations of motion are obtained: 



Nm) (^h^""^ - + (^Nm{Nm - 1)) |^^|' + Xam (iV„iV„) |0„| 



H — ( c'„ 



V2 



and 



n^'^-^^t) - e{t) ■ l] C{t) 



.ac(t) 

' dt ' 



(Al) 



(A2) 



The three time-dependent Lagrange multipliers /Xa(^), A*m(^) and e{t) appear therein. How 
to eliminate them? 
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It is straightforward to eliminate e{t). This is done by transforming the expansion coef- 
ficients as follows: 

C{t) = e-'^'^''^'^''Cit). (A3) 

Substituting Eq. (A3) into (A2) and removing at the end the 'bar' from all quantities we 
obtain: 

H^'^-"^\t)C(t)^i^^. (A4) 
at 

From Eqs. (A2-A4) we see that the role of the Lagrange multiplier e{t) is that of a (redun- 
dant) global time-dependent phase of the many-particle wavefunction ^(i) = e~*-^'^^*'^''*'^(i). 
We note that Eq. (Al) is not affected by the transformation (A3) because reduced density 
matrices are "insensitive" to a global phase of the wavefunction, namely (^(t) | . . . | ^(t)) = 

{m\---\m)- 

The next step is to eliminate the remaining Lagrange multipliers fJ^a{t) and fim{t). Making 
use of the orbitals being normalized and taking the respective scalar products of Eq. (Al) 
with {(pal and {ipm\, we obtain: 



limit) = In,, 



at 



NaNrn) {Mm \ Mm) + V2Xcon (hlh\cm^ (0a I V'm) , 

A 
2 

A 



4>r 



U"^^ - i— 



m I 0m 



> + 



+ Km {^NaNmJ {Mm \ Mm) + ) {^m I 0a> • 



(A5) 



34 



Substituting Eq. (A5) into Eq. (Al), employing the identities: 



\(pa) + V2Xcon (blblcmj 0* \^m) " ^a{t) \(f>a) = 
'a )( 0a|) I (iVa) (h^""^ - i^^ + A„ (7V„(7V„ - 1)) \<j>af + 



^m) (/i^'") - + A^ (iV^(iV„ - 1)) + 



m / lYa] 



+A ( 

+Km (NaNm) \<f>a\^ 



l^m) + (^cl^baba'j (j)a |0a) - ^m{.t) [tpm) = 

(iv^) - + A^ (iv^(iv^ - 1)) iv^r + 



(A6) 



and dividing the result, respectively, by {^^aj and ^A^y; the equations of motion (16) are 
obtained. 

To eliminate the projection operators in front of the time derivatives in Eq. (A6), i.e., 
p(«) = andP("*) = l-\ipm){'ipm\ on the left-hand sides of Eq. (16), we exploit the 

invariance property of the wavefunction ^'(i). Consider the following phase transformations 
of the orbitals and coefficients: 



0„(r, t) = e+^'^'^W0„(r, t), i^Jv, t) = e+^^-WV'„.(r, t), 
Cp{t) = e-'[(^-'^)^"W+f7^(t)]Cp(t), p = 0, . . . , [A^/2] . 



(A7) 



Combing these phase transformations, the wavefunction does not change: |^(^)) = 
T,™Cp{t)\N -2p,p-t) = T,™Cp{t)\N -2p,p;t). We should also recall that trans- 
forming the orbitals goes along with transforming the corresponding annihilation, creation 
operators (another way to look at this is that the field operators *a(r) and *m(i*) are time- 
independent, basis-set-independent quantities and, consequently, transforming the orbitals 
requires the reverse transformation of the annihilation operators). Thus, Eq. (A7) implies 
also the phase transformations 

kit) = e-^«W6„(i), tm{t) = e-^-Wc„(i) (A8) 
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for the atomic and molecular annihilation operators. 

Now, plugging Eqs. (A7-A8) into the equations of motion (16) for the orbitals and (22) 
for the expansion coefficients [see also Eqs. (A4,A6)], and choosing the phases 

Pa{t) = J i {Mt') I Mt')) dt', lm{t) = j I (ipmit') | 1pm{t')) dt' , (A9) 

equations of motion (20) and (23) are found. We note that the phases /3a{t) and ^rn{t) a-re 
real quantities since (0a(^) | 0a(^)) = 1 and {ipmit) \ i^mit)) — 1, respectively, for all times. 



2. The general multiconfigurational theory 

When the expressions for the field operators ^a(r) and 4'm(r) in terms of the time- 
dependent orbitals, Eq. (5), are substituted into the generic many-body Hamiltonian (1-3), 
one obtains: 

k,q k,s,l,q k',q' k' ,s' ,1' ,q' 

+ E <tbib,cic, +^y: i^i'kr'i'b^K + . (Aio) 

k,k',q,q' k',k,q 

The one-body, two-body and conversion matrix elements appearing in Eq. (AlO) are given 
by: 

4? = / <Pl{r,t)U'\r)<f>,{r,t)dr, 
wSji - 11 0:(r,t)0:(r',t)iy('^)(r,r')0,(r,t)0,(r',t)cirrfr', 

WS,,1' - II ^:.(r,i)V'Mr',i)iy(-)(r,r')V'.'(r,t)V^'(r',t)drdr', 
= // 0:(r,t)#(r',t)iy('^™)(r,r')0,(r,t)^,,(r',t)drdr', 
^i?"^ = //v^:,(^,i)l^(^"Hr,r')</>.(r,t)</.,(r',t)drdr' = 

= II rk'ir, i)iy(^"^"^) + " - '^'^ + ^' ^) '^'^ - ^' ^) 
= ll<Pl(v,tmr',t)W^--'^\ry)^^,,^^ (All) 
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The change of variables used for wj^J^^^^ is needed in order to perform the variation of this 
term with respect to the molecular orbitals, see below. 

Now, the expectation value appearing in the functional action (30) when expressed ex- 
plicit with respect to the orbitals reads: 



M' 



d_ 
di 



M 



(a) 



k,q=l 



(a) 



kq 



1 ^ (a) (a) 
+ 9 ^ Pkslq^ksql + 



k',q'=l 



(m) 



k'q' 



M' 



k,s,l,q=l 
M M' 



M M' 



^ Q 2-^ Pk's'l'a'^^ k's'q'l' ^ / , / , Pkk'qq'^^kk'qq' ^ 

■1 



k',s',l',q'=l 



^k'kq + Pr'-'-' ■ 



-'qkk' 



qkk' 



k,q=l k'=l 



k,q=l k',q'=l 

[N/2] 



p=0 n'P,mP 



dt 



(A12) 



where 



(a) 



dt 



dr. 



(A13) 



dt 

Equating the variation of the functional action (30) with respect to the orbitals to zero, 
making use of Eq. (A12), the following equations are obtained: 



M 

E 

q=l 



M' 

E 

q'=l 



p'i [h"'' - 4 + E ptlq^^' + E ^ 

^ ^ s,l=l k',q'=l 



+ 



M' 



k'=l 

d 



M' 



pti ( h^-^' - ) + E pK/'X^^ + E ptZK"" 



s',V=\ 



M 



M 

M 

li^q') + 

k,q=l 

M' 



■'kq 'n""=E/^Sl^/)' k'^l,...,M'. 



(A14) 



_^ J_ „(2a^'«)w(2a^rn) 

^^k,q=l j'=l 

One delicate point in performing the variation of wj^^^^ [the first term in the last line 
of the expectation value Eq. (A12)] with respect to the molecular orbitals V'^,(r,t) is worth 
mentioning. To perform this variation, a change of the integration variables r, r' to the 
center-of-mass R = and relative f = r — r' coordinates is required. Assigning thereafter 
back R — > r, f — > r', the matrix element wj^f^^ is re- written in a form, see Eq. (All), 
which is amenable to explicit variation with respect to ipl,{r,t). 

The next step it to ehminate the Lagrange multiphers {//^"^(i)} and {A*^y (^)}- Making use 
of the orthonormality properties of the atomic and molecular orbitals, {(pkit) I 0g(^)) = ^kq 
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and {ipk'{t) I i^q'it)) — ^k'q'i and taking the corresponding scalar products of Eq. (A14) with 
respect to the orbitals, we obtain exphcit expressions for the Lagrange multiphers: 



M 



(a) 



q=l 
M' 



d_ 
di 



(a) 
M' 



M 



s,l=l 



k',q'=l k'=l 



M' 



q'=l 



htj - ( t- 



J q 



d 



(m) 



M 



M' 



^ f^k's'l'q'*^^ j's'q'l' ^ 



s',l'=l 



M ^ ^ M 

(am) T^iarn) [ , J_ \^ (2a-m)T^(S 

^ Pkk'qq'^^ kj'qq' ( ^ /o i' 

fe,g=l ) k,q=l 



Substituting Eq. (A15) into Eq. (A14), making use of the identities: 

M M' 



M 
q=l 



^^"^-^^ +E^SXr^+ E pKXv 

^ s,;=l k',q'=l 
M' ^ M 



fe'=i 



u=l 



M 



M 



M' 



k'=l 



u=l / q=\ 

M' 

\- (am) u^iam] 
^ l^kk'qq'^^ k'q' 

k',q'=l 

M' 

E 

q'=l 

M 

+4E^i'"^<"^"^^-E/^K'i^"' 



|0,) + ^E/'ir.7^"^<'^'''^ 



M' 



M 



pi?' ( ^^"^^ - ) + E pK..'^^^ + E 



s',l'=l 



k,q=l 



k = l,...,M, 

\^q') + 



M' 



V2 



M' 



k,q=l 



M' 



i-Ei^"')(^-'i) 1 E 

q'=l 



u'=l 
M 



u'=l 



Pk'q' 



f)+ I Z^ Pk's'Vq'^^ s'V ^ 
s'l'=l 



+ E pK^^"^ 



kq 



k,q=l 



M 



(A15) 



A;' = 1,...,M', (A16) 



k,q=l 



and multiplying the result, respectively, by the inverse of the reduced one-body density 
matrices and summing over Yl,kLi and Yl!k'=i obtain equations 

of motion like (33) with general interactions. 
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Finally, to eliminate the projection operators P*^"^ = 1 — Ylu=i \4>u){4'u\ and P^"*^ = 
1 — X^fLi in front of the time derivatives [see Eqs. (A16) and (33)], we employ 

the invariancc properties of the wavefunction ^'(t) = \l/(t), where {0fc(r, t)} |0^(r,t)}, 
{ipk'i^^i)} — ^ {i^k'i^^t)], and {C^jp^p(t)} — > {C^^P^(^)}. For this, consider the following 
time-dependent matrices: 

DW(t), = . {Mt) I , D(™)(t), Di™) = . (^..(t) I Mt)) ■ (A17) 

The matrices D(")(t) and D("*)(t) are Hermitian matrices [because the respective orbitals are 
normalized and orthogonal to one another, {(j)k{t) \ <t>q{t)) — ^kq and {ipk'it) \ i^q'it)) — Sk'q'] 
and hence can be diagonalized 

|TW(i)}^D(")(i)T(")(i) = d(»)(i), {T^^\t)y B^^\t)T^'^\t) = d("^)(i), (A18) 

where d*^") (t) and d*^™-* (t) are the diagonal matrices of the respective eigenvalues. Now, we 
define the unitary transformations (which are symbolically integrated) : 

M 

^U^^{t)--J2D^^mi;\t) =^ U(»)(t) = e+^/'°'"'(*')^*'u(«)(0), 

k=l 
M 

iU^:rJ{t)^-Y,Df^,{t)U^;r){t) =^ U("^)(i) = e+^^*°^"'^(*')'^*'U("^)(0), (A19) 
fe'=l 

with the initial conditions defined in the hmes r — > as (see in this respect Ref. [33]): 

U(«)(^) = T('^)(0)e+*"'^''"(°), U('")(t) = T(™)(0)e+^^'*^"^(°\ (A20) 



Then, the unitary transformations of the orbitals 

M 



0,(r, t)^Yl ^kq it)Mr, t), g = 1, . . . , M, 

k=l 
M' 

V^,,(r, t) = J2 ui7;}{t)M^, t), q' = l,...,M' (A21) 

k'=l 

lead to the desired result - equations of motion (37) and (44) - where the projection operators 
P*^"-* and P^"^) have been eliminated from the left-hand sides. 

The transformation {CnPmp(^)} ^ {CnPmp{t)} accompanying Eq. (A21), carries equations 
of motion (39) for the expansion coefficients to the respective final result, Eqs. (40) and (47). 
It is instructive to obtain this result by proving that the equations of motion for the ex- 
pansion coefficients are form-invariant. Namely, if T-i^'^°'^'^\t)C{t) = are satisfied 

39 



for the untransformed quantities [{CfiPmp{t)}, {0fe(r, t)} , t)}] then T-C^^'^^"^\t)C{t) — 

i^^^ are satisfied for the transformed ones [{C^p^p(t)}, {0^(r, t)} , {■?/'^,(r, t)}] . The 
proof is straightforward. Equating the variation of the functional action (30,32) 
with respect to the expansion coefficients to zero, the result can be written as fol- 



lows: {fV^ fhP] t 
tions ^{fiP,mP;t\ 



, Vp, n^, 772.^. Since, the transformed configura- 
are given as linear combinations of the untransformed configurations 
{(n^, m^; the operator ^Pa— _ does not depend on the orbitals, and = 
we immediately get: {nP,mP-,t\H^'^°'^"^^ — |\E'(t)) , Vp, n^, m^, which concludes 
our proof. To our needs, since the transformed orbitals (A21) obey the differential 
conditions I — 0;^;? — 1, • • • , M and (jPk'l'^q''^ = 0;k',q' = 1,...,M' [see 
Eq. (38)], the respective equations of motion for the transformed coefficients boil down 
to H^'""'"'V)C(t) = [see Eqs. (40) and (47)]. 

APPENDIX B: MATRIX ELEMENTS WITH MULTICONFIGURATIONAL 
WAVEFUNCTIONS IN SYSTEMS WITH PARTICLE CONVERSION 

There are two types of matrix elements in the theory. The ffist type are matrix elements 
of the many-body Hamiltonian with respect to the configurations. These matrix elements 
are expressed using the matrix elements of the one-body terms, two-body interaction terms, 
and the conversion term with respect to the atomic and molecular orbitals. The second type 
of matrix elements are the matrix elements of the reduced density matrices appearing in the 
theory, which are expressed in terms of the expansion coefficients. 

In this appendix we prescribe these matrix elements. It is easy to connect the matrix 
elements of particle-conserving operators to the corresponding matrix elements appearing 
in the available multiconfigurational theories for identical particles and mixtures. This 
assignment will shorten substantially the discussion below. The matrix elements of particle 
non-conserving operators are new and will be presented in full details. 

1. Matrix elements of the Hamiltonian 

The many-body Hamiltonian (1-3) is written as a sum of particle-conserving and particle 
non-conserving parts: — ^(""i) -\- |y(2a^TO) xhe matrix elements of the particle- 
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conserving part H^"-"^\ see Eqs. (1,2), between two general configurations derive from the 
following relation: 



fF, wF\ t 



d 
di 



'P'.- 



jp,p' 



fF, rrV; t 



d_ 
di 



'p. 



(Bl) 



Thus, it corresponds to a matrix element of a mixture with N — 2p atoms and p molecules 

without conversion. The matrix elements of the Hamiltonian of a mixture of two kinds of 
bosons with respect to two general configurations have been prescribed within the respective 
particle-conserving multiconfigurational theory for Bose-Bose mixtures, the MCTDH-BB 
theory, see Ref. [54]. 

To evaluate the matrix elements of the particle non-conserving part of the Hamiltonian, 
W^'^"-^'^'^ = j^]Y{rn-.2a) ^ (1^3)^ between two general configurations we have 

to introduce a shorthand notation for different configurations. The reference configura- 



tion is denoted by |n^,m^;t) 



1) • • • I ""-fci • • • ) 



p 



: m 



15 



,mi,,...,m' 



%-t). We 



remind that the occupation numbers satisfy the relations: |n^| — n\ + . . . + n\j — N — 2p 
and |m^| = -|- . . . -|- m^, = p, where p is the number of molecules. Now, the config- 
uration |fi^~\ m^r^ t) = \rfi, ... ,nl + 1, ... ,nP + 1, ... ,711^ : rrfi, ... ,ml, - 1, ... , m%] t) 
differs from |n^, wP; t) by having p — 1 molecules and N — 2p -\- 2 atoms, where 
a molecule in the A;'-th orbital has dissociated to two atoms, one in the A;-th or- 



bital and the second in the g-th orbital; and the configuration 



\n' 



p-i 

kk 



, m 



1^1' ■ ■ ■ , + 2, . . . , : 



m 



1' 



, m^, — 1, ... , m^,; t) differs from IfP, mF\ t) by having p — 1 



molecules and N — 2p + 2 atoms, where a molecule in the k'-th orbital has dissociated to 
two atoms, both in the k-th orbital. We employ a nomenclature where the same ordering of 
the orbitals (pi, ... , (pM and ipi, . . . , ipM' as in Eq. (26) is kept for all configurations. In this 

nomenclature the following states arc equivalent: \ni^'^,'i7i'l7^',t^ = \n^k^ , iTi^j''^ ', t^ 

With this notation, the non-vanishing matrix elements of the particle non- conserving part 
of the Hamiltonian follow from 



(^fF, wF; t 



fF, nF; t 



|y(2o^m) 



p—1 -fp—l_ 



''^kq ' ''^k' 
->p—l 



^^k'kq 



1 

7^' 



' / — fT^^^k'kk 



(2a^m) 



ml,{nl + l){nl+l) , k < q, 
ml,inl + l)inl + 2) , (B2) 



and the relation ('nP,rfF;t 



fF,ffF;t'^Y ■ 



We note that W^Sg""""^ = ^Sq^"^^ because W^(2«--"^)(r, r') = W^'^^^'^\y' ,y). To summa- 
rize, direct coupling in the matrix representation of the Hamiltonian with respect to the 
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configurations exists due to the particle non-conserving part of the Hamiltonian 
between configurations with p and p — 1 molecules only, for p = 1, ■ ■ ■ , 



2. Matrix elements of reduced density matrices 



The multiconfigurational ansatz (26) can be written in the following form, 

[N/2] 
p=0 



(B3) 



where each |\l/p(t)) = 'Yli-nv mp^nvmp{t) \TP'-,'fn^\t) is a (non-normalized) many-particle wave- 
function with a definite number of p molecules and N — 2p atoms, and thus "describes" 
a bosonic mixture without conversion. Consequently, the matrix elements of the particle- 
conserving reduced density matrices can be expressed as follows: 

[N/2] 

p=0 

[N/2] 



cliCq' 



p=0 

[N/2] 

Ait) = {bibikb,) = J2 , 

[N/2] 

(m) 



p=0 

[N/2] 

p=0 



blbgcl,Cg< 



(B4) 



In other words, the matrix elements of the particle-conserving reduced density matrices 
can be readily read from the reduced density matrices of the respective particle-conserving 
multiconfigurational theory for Bose-Bose mixtures, the MCTDH-BB theory, see Ref. [54]. 

The matrix elements of the particle non-conserving reduced density matrices are given 
explicitly by 

[N/2] 



(2a^m) 
Pk'kq 



{2a^m) 
Pk'kk 



clbkbq) = E E ^nP,m^C!^l-\rnl7' \/ K'K + l)« + l) ' ^ < 

p=0 nP,mP 

[N/2] 



p=Q fiP,mP 
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All other matrix elements are derived from the symmetry of the conversion operator, 
Pu'lk"^^ = pSr^"^' Hermiticity relation pj^fe7^"^(t) = ■ 
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